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ABSTRACT 

Successful IUE observations of the equatorial sunlit atmosphere of Jupiter and Saturn have been 
obtained. Spectra containing atomic and molecular hydrogen and solar reflection continuum emis- 
sions have been analyzed, with the purpose of determining the long term temporal behavior of the 
electroglow process. Quantitative estimates have been established for the first time using a model 
analysis of the short wavelength region of the spectrum. Both systems show varying degrees of long 
term variability in hydrogen emission rate, but the time scale is too short to determine whether there 
is a dependence on solar cycle activity. As part of the emission modeling program, a preliminary 
point source spreading function for the IUE SWP instrument has been established, suggesting a 
wavelength dependence in spectral line width different from previous analyses. Further IUE obser- 
vations are planned for both Jupiter and Saturn. 

INTRODUCTION; OBJECTIVES 


JUPITER 

A question of recognized importance to the understanding of the source and flow of energy in 
Jupiter’s atmosphere, is the long term behavior of the EUY emissions. The initial epoch for usable 
observations relevant to this question is ~ 1 970. The principle reference points are rocket observa- 
tions in 1968-72, Pioneer 10 (P10) and Voyager (V) spacecraft encounters in 1973 and 1979 near the 
times of solar minimum and maximum activity respectively, and IUE measurements beginning in 
1978. Recent reanalysis of rocket observations obtained originally by Giles et al. (1976), together 
with the P10, Voyager, and IUE data indicate that the H 2 band emissions have remained relatively 
constant, while the H Lyot emission has varied in intensity by an order of magnitude (Shemansky 
and Judge 1986; Appendix A). The more recent IUE measurements in this series (1985-1986) which 
have now provided a measure of both H Lya and H 2 emissions, with the application of spectral 
models, suggest that the H 2 bands may have increased in brightness at solar minimum whereas the H 
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Lya values decreased. This result contrasts with H Lya measurements of Uranus over the 
1982-1985 period in which there is no long term variation (Clarke et al. 1986). Table 1 shows a 
summary of results. On this basis Jupiter exhibits a substantially larger variation at the H Lya 
emission rate than the solar source between minimum and maximum activity. The H Lya bulge 
phenomenon showing longitudinal variation in emission rate is unique to Jupiter. The proposed 
bulge mechanism (Shemansky 1985) requires particle energy deposition in the exosphere, a symptom 
of which is the observed equatorial H 2 Rydberg bands. The apparent strong variability of the H 
Lya emission both in absolute terms and in relation to the H 2 Rydberg systems, is basically not 
understood, but it appears to imply a changing abundance of H forced primarily by variations in 
loss rate rather than production rate. 

The April 1985 and November 1986 IUE results obtained in this program show a reduced abso- 
lute H Lya intensity, but not to the extent obtained from the equivalent results near solar minimum 
in 1974. The analysis of the 1985, 1986 measurements also have obtained a measure of H 2 Rydberg 
brightness which may have increased relative to 1979. The reason for an apparent decreased abun- 
dance of atomic hydrogen in the face of an increase in the major dissociation source is an interest- 
ing mystery. Further investigation is required to determine the extent of the dependence on the 
solar cycle. 

It should be emphasized that although the brightness of the strong Jovian H Lya line has been 
monitored since shortly after the launch of IUE, little of this data can be used to establish the cor- 
responding intensities of the H 2 band emissions. These emissions are complementary to the H Lya 
intensity, providing a spectroscopic diagnostic of the dissociation and energy deposition rates. 
Understanding the role played by solar EUV radiation in catalyzing this emission will require 
knowledge of the behavior of the Lya and H 2 band components throughout the solar cycle. 


SATURN 

We propose to obtain further observations of Saturn’s sunlit atmospheric emissions, on the basis 
of preliminary analysis of very recent measurements obtained in the 9th year IUE program on 1986 
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DOY 252 and 254. These results show a mean H Lya brightness a factor of 5 below the value 
measured by Voyager 1 in November 1980. The details of the analysis indicate that the IUE and 
Voyager results can be directly compared as discussed below, and we conclude that the apparent 
change in emission brightness is real. The present result is remarkable because the disk averaged 
brightness of Saturn in H Lya (0.88 +_ 0.2 kR) is less than the mean value obtained for Uranus by 
Clarke et al. (1986). The Saturn disk averaged brightness in H Lya radiation has not been moni- 
tored on a frequent basis so that the time scale of the variation is not well defined. Since the 
launch of the IUE observatory and the Voyager encounters with the outer planets it has become 
apparent that, contrary to our expectations, the solar UV flux contributes only a small fraction of 
the energy deposited in the upper atmospheres of Saturn and Jupiter (Broadfoot, et al, 1981a, Broad- 
foot et al., 1981b; Shemansky, 1985; Yelle, et al., 1986. Evidence for this is the high exospheric 
temperatures (1000 K for Jupiter, 420 K for Saturn and 800 K for Uranus), the intense, non- auro- 
ral, H 2 band emissions (3 kR from Jupiter and 1 kR from Saturn) and the H Lya to H 2 band 

brightness ratio which indicates a high altitude, collisional source for both emissions. These devel- 
opments imply that the measured solar UV flux can no longer be used as the sole input for calcula- 
tions of the physical state of Saturnian and Jovian upper atmospheres. Rather, we must rely on the 
H 2 band brightness to measure the energy deposition and infer the exospheric temperature, ioniza- 
tion rate and atomic H mixing ratio in the upper atmospheres (Table 2 shows estimates of dissocia- 
tion rates based on Voyager data). Fortunately our understanding of the properties of the H 2 
molecule, and consequently synthetic spectra, have advanced sufficiently to allow the accurate deter- 
mination of total emission brightness, energy deposition and electron temperature (given a sufficient 
signal to noise ratio) from the measured H 2 band spectrum (Shemansky et al., 1985; Shemansky, 
1985). Similarly, the H Lya line is produced predominantly by collisional excitation rather than 
resonant scattering of solar Lya. The resonant scattered contribution is less than 50% on both 

Saturn and Jupiter [Shemansky, 1985; Yelle et al., 1986]. For this reason the H Lya to H 2 band 

brightness ratio directly reflects the atomic H mixing ratio and hence the altitude of the excitation 
source. The H Lya to H 2 band brightness ratio provides an important clue to the nature of the 
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underlying energy source and the physical state of the upper atmosphere. 

The relative instability of Jupiter and Saturn in apparent abundance of atomic hydrogen must 
carry clues to the source process if the nature of the variability can be identified. It is not clear 
what the relationship to solar cycle is, because all three outer planets, Jupiter, Saturn, and Uranus 
according to available data show different behaviors in relation to solar activity. 

RESULTS; DISCUSSION 


JUPITER 

We now have an IUE point of reference in August 1984, April 1985, and November 1986. The 
analysis of this data (Table 1, Figure 1) shows a moderately reduced H Lya intensity, but a relative 
bulge amplitude roughly the same as that observed by Voyager. This appears to imply that the H 
Lya is excited at about the same atmospheric and ionospheric density levels as in 1979. The solar 
reflection continuum is at the same level as that obtained by both Voyager and IUE observations in 
1979 with the exception of the 1986 result which shows a value 16% larger (Table 1). H 2 band in- 
tensities obtained from model calculations compared to the April 1985 IUE data indicate values that 
are basically the same as an analysis of the Clarke, Moos and Feldman (1982) results, I s (H 2 
Lya+Wj^) = 4.4+J.6 kR. This value is ~45% larger than the Voyager measurement in 1979. Ana- 
lysis of the November 1986 spectrum (Figure 2) gives a value 5.4 +_1.7 kR, moderately higher, than 
the 1985 and earlier value, but about a factor of- 2 larger than the Voyager value of I s (H 2 
Lya+Wft) = 3.0 kR. We have proposed further IUE observations of Jupiter to extend the data base 
further into the solar minimum period. The additional data we are proposing to obtain here would 
be composed of a series of bracketed SWP exposures designed to measure H Lya as well as the 
fainter H 2 features. These data would characterize the Jovian equatorial emission under solar mini- 
mum conditions. The importance of obtaining an equatorial H 2 emission spectrum lies in the fact 
that the emission is clearly unaffected by hydrocarbon absorption, and provides a direct measure of 
energy deposition rates. 
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SATURN 

We have obtained measurements of the Saturn electroglow spectrum on 1986 DOY 252 and 254 
in this program. Observing conditions were remarkably good. Analysis of the data requires rela- 
tively more model analysis than is the case for Jupiter, in order to account for extinction effects, 
geocorona, and non uniformity of the source in the field of view. 

Table 3 shows a comparison of earlier results. The corrections for foreground and background 
LISM contributions and foreground geocorona require model calculations that have not yet been 
made for the Clarke et al.(1981) observations. Moreover at the time of these observations, Saturn 
was not in a favorable position in respect to extinction effects in the intervening LISM. Crudely 
estimated extinction coefficients are given in Table 3. The measured H Lya brightness from IUE is 
substantially lower than the VI result, which we attribute mostly to extinction, by the interplanetary 
medium (ISM/IPM) partly to the analysis method, and possibly partly to variability in the source. 
The 1980 IUE observations are close enough to solar maximum to constitute a solar maximum refer- 
ence. Observations in 1986 and 1987 will provide measurements near solar minimum, under much 
more favorable conditions (see Table 3 for extinction estimate in 1986). 

It should be emphasized that although the brightness of the H Lya line on Saturn has been 
monitored over much of the lifetime of IUE, little of this data can be used to establish the corres- 
ponding intensities of the H 2 band emissions. As we have pointed out H 2 emissions are comple- 
mentary to the intensity of the H Lya line, providing a spectroscopic diagnostic of the excitation 
process. Understanding the role played by solar EUV radiation in stimulating this Saturnian EUV 
emission will require knowledge of the behavior of both components at conditions of solar maxi- 
mum and minimum. Unfortunately the observations designed to provide H 2 band intensities in this 
program have not been reduced to this extent for lack of available time (1986 DOY 252, 254) for 
analysis. We therefore do not have a measure of how accurately the band emissions can be meas- 
ured. 

The observational sequence in the last program (1986 DOY 252, 254) is shown in Table 4. The 
Table includes the reduction of the H Lya brightness. We regard this as the best combination of 
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exposures required to obtain both H 2 band H Lya and solar reflection continuum. The reflection 
continuum albedo is low enough that moderate changes in eddy diffusion coefficient may be detect- 
able. 

We require models of the LISM as well as the geocorona in order to properly reduce the data. 
We have developed a model of the geocorona for this purpose, and we have used the recent Ajello, 
et al. (1987) model of the LISM to take complete account of the extraneous H Lya signal in the 
data. The result of the application of the models is shown in Table 4. This is the only IUE obser- 
vational data from Saturn that has been analyzed in this way, and it will be necessary for compara- 
tive purposes to reanalyze the earlier data obtained by Clarke et al. (1981). Table 4 shows the 
sequence of observations with the measured H Lya intensity in Col. 4. Interspersed with the Saturn 
observations are background data obtained 60" north of the planet. The background is composed of 
geocoronal and LISM components as estimated in Col. 5 and 6 of Table 4. The component of the 
LISM in the foreground to the planet is estimated to be 0.463 kR as given in Table 3, using the 
Ajello, et al. (1987) LISM model. The total LISM intensity in the look direction was 0.622 kR. 
Although the observed emission brightness varied from -3.1 kR to -1.8 kR during the sequence, 
most of the variation appears to be caused by the geocoronal component, and the derived Saturn 
emission brightness is basically constant during the sequence with an estimated mean value of IHLya 
= 0.88+_.2 kR. The geometry for the 1986 DOY 252, 254 observations is ideal in the sense that the 
planet is essentially directly upstream from the earth relative to the bulk inflow of the LISM, so that 
extinction of the planetary signal is at a minimum. However, the observing geometry at the time of 
the Clarke et al. (1981) observations indicates a substantially larger extinction factor (Table 3). If 
these estimates are correct the evidence indicates that the Saturn emission brightness in H Lya has 
varied by a factor of 5. The inference is that the abundance of atomic hydrogen on the planet must 
have changed significantly. 

The analysis of the 1986 DOY 252, 254 data is incomplete, because there has been insufficient 
time to reduce line by line data to obtain estimates of the H 2 equatorial emission. We therefore do 
not have a good assessment of the uncertainties in determining H 2 band intensities. In contrast to 
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the H Lya analysis described above, auroral and ring reflection components must be removed from 
the data in order to analyse the H 2 band electroglow emissions. 


DEVELOPMENT OF EMISSION MODELS 


We use an H 2 emission model that has been employed extensively in the analysis of both Voy- 
ager observations of H 2 emission from Jupiter and Saturn (Shemansky and Ajello, 1983; Shemansky 
1985) as well as interpretation of electron excited H 2 emissions in the laboratory (Ajello et ak, 1984, 
Shemansky et ah 1985). We have developed and applied an accurate model of the IUE instrument 
transmission function for use with our calculations of the emission structure of the H 2 Rydberg sys- 
tems (see accompanying report). It is composed of a convolution of a gaussian (point source) and a 
trapezoidal function. The FWHM of the trapezoidal function is constant as a function of wave- 
length. However on the basis of our analysis of a point source H Lya line, and results reported by 
Cassatella et al. (1983), we apply a gaussian function with a linear FWHM dependence on wave- 
length, ranging from 3.44A at 1216A to 6.40 A at 1900 A. This analysis differs from that of Cassa- 
tella et al. at short wavelengths (see accompanying report). The instrument transmission function is 
clearly an important factor in the analysis of discrete structure, and the variation of the gaussian 
point source function measurably affects results. We use an H 2 emission model that has been emp- 
loyed extensively in the analysis of both Voyager observations of H 2 emission from Jupiter and 
Saturn (Shemansky and Ajello, 1983; Shemansky 1985) as well as the interpretation of electron 
excited H 2 emissions in the laboratory (Shemansky et ak, 1985). As far as we are aware it is the 
most comprehensive model of H 2 emission employed so far. 

Figure 2 shows the analysis of the November 1986 spectrum. All of the spectral components 
are modeled, including solar reflection continuum. 
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TABLE 1 


Jupiter sub-solar brightness with no correction for 
ISM/IPM extinction 


Obs. 

Date 

I<H $o re) 

V 

I(H Lyct) 

I (H-> Ly+WR) 

(1600A) 

R/A 

III 

CML 

Remarks 

Rocket a 

1972/144 

2.4 

1 . 

6.4 

50. 

100 

Disk Average 



3.1 

1.5 

0.5 

74. 

100 

Calculated 








Sub-solar Point 

Rocket & 

1978/335 


13 



106 

Disk Average 


1978/335 


19 



106 

Sub-Solar Point 

IUE c 

1979/158 




74. 


Sub-Solar Point 

V2 d 

1979/187 

3.0 

16 

5.3 

76. 

240-330 

Sub-Solar Point 

V2 d 

1979/187 

3.0 

22. 

7.2 

65. 

60-150 

Sub-Solar Point 

IUE « 

1979/120 


15 



100 

Sub-Solar Point 


-150 







IUE e 

1980/120 


12 



100 

Sub-Solar Point 


-180 







IUE f 

1985/99 

4.4+_l .6 

8.0 

1.8 

76+6 

100 

Sub-Solar Point 

IUE f 

1985/99 

4.4+_1.6 

5.8 

1.3 

76+6 

280 

Sub-Solar Point 

IUE f 

1986/304 


8.2 

1.6 


100 

Sub-Solar Point 

IUE f 

1986/304 

5.2+.1.5 

6.0 

1.1 

85+_6 

280 

Sub-Solar Point 


a) Judge and Shemansky (1985) analysis of Giles et al. (1976) rocket experiment, 

E-W H bands distribution shows no limb darkening (Shemansky, 1985) 

b) Clarke et al. (1980) rocket experiment 

c) Present work, estimate from SWP 5448 

d) Shemansky (1985) analysis of V2 data 

note: I(1600A) - 0.43; Table 3 should read "Atmospheric Reflection 1660 A; Solar" 
1(1 660 A) 

e) Clarke et al. (1981a) corrected upward by a factor of 1.16 on 1983 correction 

to IUE aperture size. 

f) Present work Holberg and Shemansky (1985). 




TABLE 2 

Exobase Production Rates and Escape of Atomic 
Hydrogen from the Outer Planets 2 



Jupiter 

Saturn 

Uranus 

(1) 

1.6 x 10 30 

3 x 10 29 

2 x 10 28 

(2) 

5.1 x 10 9 

1.7 x 10 9 

5 x 10 8 

(3) 

19. 

6. 

2.3 

(4) 

0 

~10 

i 

OO 

o 

(5) 


1.6 x 10 9 

8 x 10 8 

(6) 


3 x 10 28 

1.6 x 10 28 

(7) 


10 11 


(8) 


2 x 10 4 

2 x 10 4 

(9) 

5-7.3 

2.8-5.3 

>2. 

(10) 

3.0 x 10 9 

1.0 x 10 9 

4.0 x 10 8 


a These estimates exclude auroral activity (see text). 

(1) Dissociation production rate (s -1 ); 

(2) (cm2 s -i) 

(3) HI escape energy (eV); 

(4) Escape yield (%) 

(5) Loss rate (kg yr -1 ); 

(6) (atoms/s) 

(7) Loss lifetime for 200 km am H 2 (yr) 

(8) Exosbase ion density (cm -3 ); 

(9) Y IS (Eq. 1) 

(10) Ionization Rates (ions cm -2 s -1 ) 
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TABLE 3 

Rough estimate of extinction by the 
ISM/IPM for observations of H Lya 



a 

Date I(H Lya) 
(kR) 

b 

r 

c 

Ic(H Lya) 
(kR) 

Ic(H Lya) 
I(H Ly+WR) 

A III 




Jupiter 




Rocket 

1972/144 1.5 

0.236 

1.9 

0.61 

100 

Rocket 

1978/355 19. 

0.30 

26 


100 

V2 

1979/187 22. 

0.0 

22. 

7.2 

60-150 

IUE 

1980/120 15. 

0.30 

20. 


100 


/ 1 50 





IUE 

1980/120 12. 

0.30 

16. 


100 


/150 





IUE 

1985/99 8. 

0.236 

10. 

2.3 

100 


IUE d 

1980/19 

0.9 

Saturn 

.85 2.2 


IUE d 

1980/125 

0.7 

.85 

1.5 


VI 

1980/316 

4.9 

0.0 

4.9 

5.3 

V2 

1981/236 

3.0 

0.0 

3.0 

2.8 


a) Brightness before correction for extinction 

b) Calculated optical thickness (see text) 

c) Brightness after correction for extinction 

d) Clarke et al. (1981b) corrected upward by factor of 1.16 on basis 

of 1983 correction to IUE aperture size. Saturn is known to 
have no limb darkening effect in H Lya and therefore no limb 
darkening corrections are applied to the measured data. 
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TABLE 4 


1986 DOY 252 and 254 SATURN IUE OBSERVATIONS OF H LYa 


SWP 

Exposure 

Duration 

(Min) 

Scattering 

Background 

(FN/S) 

ILy-a 

Observed 

(kR) 

ILy-a 

Geocorona 

(kR) 

ILy-a 

ISM 

(kR) 

ILy-ab 
Saturn Emission 
(kR) 

29 1 70 a 

30 



0.917 

0.295 

0.622 


29172 

80 

2.83 

2.862 

1.651 

0.508 

0.976 

29173 

30 

2.65 

3.137 

2.019 

0.508 

0.849 

29 I 74 a 

15 

— 

2.276 

1.593 

0.622 

— 

29175 

15 

2.25 

2.295 

1.220 

0.508 

0.788 

29176 

25 

2.83 

1.979 

0.920 

0.508 

0.765 

29177 

25 

2.50 

1.747 

0.686 

0.508 

0.768 

29 1 89 a 

30 

— 

2.166 

1.604 

0.622 

— 

29190 

30 

2.39 

2.220 

1.124 

0.508 

0.817 

29191 

30 

2.27 

1.903 

0.803 

0.508 

0.822 

29192 

15 

2.37 

1.819 

0.625 

0.508 

0.953 


WEIGHTED MEAN 0.88+0.2 


a - Sky Background Observation (60" N of Saturn) 
b - Implied Saturn Emission = (I 0BS - I ISM - IgeOCORNa)/- 72 
c - I(Ly-a) ISM = [I(Lya) ISMoo ]*.28 + [I(Lya) ISMSAT ]*0.72 
I(Lya)is Moo = 0.622 kR 
I(Lya)js M s AT =0.463 kR 


[ (v/hJ) I ]°'Boi 
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CML (degrees) 


Fig. 1 H Ly a brightness Jupiter 
center, calculated from 
IUE measurements in 1984- 
1986, as a function of CML. 
Geocoronal and LISM 
components have been re- 
moved using model calcula- 
tions referenced to 
interspersed IUE background 
measurements. 

PI August 1984 

© April 1985 

A November 1986 


Figure 1 
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Fig. 2 Reduced composite spectrum 

of Jupiter obtained November 
1986. The data is a combin- 
ation of SWP29576 (15 min.) 
SWP29576 (80 min.) and 
SWP29577 (90 min.) such that 
the entire spectrum contains 
unsaturated signal. The 
heavy overplotted line is a 
model calculation of electron 
excited ”2 > combined with 
an H Lya line and a model 
of the solar reflection 
continuum. 


Figure 2 
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Abstract 


Analysis of the Pioneer 10 and rocket observations of disk averaged 
emission from the sunlit atmosphere of Jupiter indicates that the spectrally 
integrated EUV brightness was reduced by at least a factor of 2 relative to 
Voyager spacecraft observations in 1979. Most of the variation is caused by 
the H Ly a component in the spectrum, which was reduced ~ 1 order of 

magnitude near the time of solar minimum in 1972-1973. Although the analysis 
of the data does not produce entirely consistent results, the weight of 

evidence points to a factor of order - 2 lower abundance of HI in Jupiter's 
atmosphere in 1972-1973 relative to 1979. The low emission rate in H Ly a 
near the time of solar minimum in this proposed scenario is caused by an 
electroglow energy deposition rate reduced by a factor of ~ 3. The apparent 
reduced abundance of HI implies a reduced thermospheric temperature, even 
under the assumption of a constant electroglow deposition rate. 


IUE SWP Point Source and Filled Field 


Spectral Transmission Function 

D. E. Shemansky, J. B. Holberg, and D. T. Hall 
Lunar and Planetary Laboratory 
University of Arizona, Tucson, A Z 85721 

In the analysis of discrete emission spectra obtained with IUE, we have found 
that it is important to apply an accurate spectral transmission function. In the 
course of applying model calculations to emission spectra of the outer planets 
using the SWP camera we have established what we believe to be a more accurate 
function particularly at the shortest wavelengths. Cassatella et al. (1985) have 
provided a valuable analysis of the system in both the dispersive and non 
dispersive coordinates. Figure 6 of the Cassatella et al. (1985) paper shows a 
plot of the FWHM for the SWP based on point source observations in the large and 
small aperatures, showing variation in the FWHM width of emission lines as a 
function of wavelength, using extracted spectra as given by the IUESIPS. 
Cassatella et al. suggested that resolution was at a peak in the 1300-1400 A 
region, although uncertainty in the measured widths increased at the shorter 
wavelengths due to the poor quality of the available sources. Many of the 
features of major interest in outer planet atmospheric observations occur in the 
vicinity of the H Lya line. In order to produce a better measure of the point 
source spreading function, we obtained a measurement of the H Lya line in emission 
from the geocoronal background, using the small aperture. This spectrum in flux 
numbers is shown in Figure 1. Although the source fills the field of the small 
aperture, the line is assumed to be equivalent to a point source function, because 
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the ultimate instrumental resolution is the limiting factor at the solid angle 

defined by the aperture* A model calculation using a Gaussian function with FWHM 
o 

=* 3.44 A provides an optimum fit to the line. Figure 1 shows the Gaussian curve 

integrated over the mean channel interval of the instrument in comparison to the 

observed line. The fit clearly appears to be exceptionally good. If one applies 

o 

the data plotted by Cassatella et al. a value FWHM - 5.3 A would be obtained, 
approximately 50% larger than the present result. On this basis with the trend 
shown in the Cassatella et al. data we recommend that a linear function be 
applied to the calculation of the FWHM values of the Gaussian shape as a function 
of wavelength, 

FWHM = a X + b , (A) (1) 

1150 A iX 2000 A, 

where 

a = 4.444 x 10 , 
b = -1.984. 

The quantity FWHM varies from 3.44 A ( @1216A) to 6.9 A (@2000A). Figure 2 

shows a plot of the proposed linear FWHM function against the data provided by 

o 

Cassatella et al., and the present result at 1216 A. In model calculations for 
spectra of sources filling the SWP large aperture, we have applied a convolution 
of the point source function given above with a constant trapezoidal shape 
function having a full width at the peak of 

FWP = 7.24 (A), (2) 

with a value 

FWHM = 9.96 (A) (3) 

The parameters for the trapezoidal function have been designed to provide a best 
fit to the H Lya line of uniform intensity filling the large aperture. Figure 3 
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o o 

shows a direct comparison, of 1216 A and 1800 A monochromatic lines modeled for a 

uniform spatially diffuse source, illustrating the effect of the variational point 

source function. The spectra in this case are given in absolute units, using the 

IUE photometric calibration curve as it would be applied to real data. Each line 

is fixed at a brightness of 1R. The Figure 3 plot of differential (R/X) brightness 

indicates differences in peak value of the order of 10%, the 1216X line being 

significantly narrower. Another effect that should be noted in Figure 3 is that 
o 

the 1216A line shows a noticeable asymmetry in its shape. This is caused by the 

fact that the IUE calibration curve is not flat in the 1150-1300 X region of the 

spectrum; the application of the calibration curve to the instrumental 

transmission function distorts the shape of the resultant calibrated spectrum. At 

12161 this causes very little error in the integrated line intensity (< 1%), but 

model analysis should follow the procedure described here in order to accurately 

o 

describe the calibrated spectrum at and below 1216 A* The FWHM value of the 1216 
X feature in this calculation is 10.5 A. Note that the calibration curve distorts 
the shape of the 1216 X line in Figure 3 sufficiently to cause the peak in the line 
profile to occur ~ 1 X shortward of the position of the source line. Attempts to 
adjust the wavelength scale using calibrated spectra may therefore introduce 
errors; such adjustments should be made to the flux number spectra prior to 
application of the calibration curve. From this point of view spectral analysis 
with or without the use of models should be done using flux number spectra, with 
a subsequent conversion to absolute quantities; in effect, the calibration process 
applies conversion factors to photons that do not match the wavelength of the 
conversion factor, causing distortion of the true instrumental transmission 


function. 
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We recommend that model calculations be carried out in a two stage process 
in which point source gaussian lines are first placed in a buffer with integration 
channel widths of 0.1 A. The buffered data is then treated as a monochromatic 
line source for the trapezoidal transmission function. The resultant convoluted 
spectrum can then be integrated to match the mean channel width of the TUESIPS 
data. A code for a normalized trapezoidal function integrated over user selected 
channel widths is given below. 


Reference 

Cassatella, A.^J. Barbero, D. Benvenuti, On the international ultraviolet explorer 
(TUE) point spread function at low resolution, IUE Newsletter #24, June 1984, p. 
84. 
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SUBROUTINE TFUN2(W.WZ.SN1, SN2.ffRES.IERR) 

DES 1/4/83 

TRAPEZOIDAL TRANSMISSION FUNCTION 
W= WAVELENGTH OF DELTA FUNCTION LINE 

WZ=WAVELENGTH OF SHORT WAVELENGTH EDGE OF 1 ST CHANNEL 
OF THE TRANSFORM DATA FILE 

SN1=FWHM/WRES , WHERE FWHM IS FWHM OF TRANSMISSION FUNCTION 
SN2=HWP/WR£S , WHERE HWP IS THE HALF WIDTH OF THE TRAPEZOID PEAK 
WRES* WIDTH OF DATA CHANNEL IN UNITS OF W 
IERR=-1 INDICATES OVERRANGE OF TRANS DIMENSION IN COMMON 
THE DATA DEPOSITED IN THE TRANSFORMED FILE IS I*TRANS 
WHERE I IS THE LINE INTENSITY 
TRANS(25) ON EXIT CONTAINS LINE TRANSMISSION FUNCTION 
IN LOCATIONS TRANS(l) THRU TRANS(NMAI) 

MINCH=FIRST LOCATION IN DATA FILE FOR SUMMATION OF THE 
TRANSFORMED LINE; SUM I*TRANS(1) INTO MINCH 
MAICH=LAST LOCATION IN DATA FILE FOR SUMMATION OF THE 

TRANSFORMED LINE; SUM I«TRANS(NMAT) INTO MATCH 
FIRST CHANNEL IN DATA FILE IS NUMBERED +1 
WAVELENGTHS INCREASE WITH INCREASING CHANNEL NUMBER 
THE TRANSMISSION FUNCTION IS NORMALIZED TO 1.00 
SN2=0.0 DEGENERATES THE FUNCTION TO A TRIANGLE 
DIMENSION UF(4) ,RF(4) 

COMMON TRANS(25),NMAX, MINCH, MATCH, NCH 
DATA RF/ 1.0, -1.0, -1.0, 1.0/ 


IERR=0 

SNCH=(W-WZ) /WRES 
NCH=TFIX(SNCH) 

PCH= SN1- SN2+NCH- SNCH 

NP=IFIT(PCH)+1 

DEL=PCH-NP 

NM=IFIT(2. 0 *( SN1- SN2 ) -DEL) 

NMAT=NM+1 

IF(NMAT.GT.25) GO TO 70 
MINCH=NCH+1-NP 
MATCH® MINCH+NM 
DO 10 1=1,25 
10 TRANS(I)=0.0 
UF(1)=0.0 
UF(2)=SN1-2.0*SN2 
UF(3)=SN1 

UF(4)=2.0*(SN1-SN2) 

TN=SN1*UF(2) 

TN=0.5/TN 
DO 60 I=1,NMAT 
N=I-1 

DO 50 1=1,4 
ALPH=N+DEL-UF(K) 

IF(ALPH .GE. 0.0) SL=2.0*ALPH+1.0 

IF(ALPH .LT. 0.0 .AND. ALPH .GT. -1.0) SL=(ALPH+1.0)**2 
IF(ALPH .LE. -1.0) SL=0.0 
50 TRANS(I)=TRANS(I)+RFtt)«SL 
60 TRANS(I)=TN*TRANS(I) 

65 CONTINUE 
RETURN 
70 IERR=-1 
GO TO 65 
END 



Figures 


Fig. 1 A comparison of the observed and modeled profile of the L y a line in 
the SWP small aperture. The observed profile represents the net 
flux number spectrum from a 335 m exposure (SWP25630) of the 
geocoronal Lya background. 

Observed data; oooo. 

Model calculation; . 

An optimal fit to the observed data is shown calculated using a 
Gaussian profile with FWHM = 3.44 A integrated over channel widths 
corresponding to the IUE data interval. 

Fig. 2 Point Source Gaussian FWHM as a function of wavelength, recommended 
for use in the analysis of IUE SWP data. The formula for the linear 
function is given in the text. 

□ - FWHM determined from Figure 1. 

O- Small aperture data from Cassatella et al.; 2-3 measurements 

• - Small aperture data from Cassatella et al.; 4-9 measurements 

Fig. 3 Model calculations of IUE SWP monochromatic lines at 1216A and 
1800A, from a spatially diffuse source. The data is calibrated (R/A) 
using the IUE SIPS calibration function. Each line has an integrated 
brightness of 1 R. The line shapes are calculated using the 

convolution of Gaussian point source and Trapezoidal aperture 

function with parameters given by Eqs. (1) - (3). Distortion and 
shifting of the peak position of the 1216A line is caused by the 
strong wavelength dependence of the calibration curve in that region 
of the spectrum. 
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